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We explore the electromagnetic gauge invariance of the hadron tensor of the Drell-Yan process 
with one transversely polarized hadron. The special role is played by the contour gauge for gluon 
fields. The prescription for the gluonic pole in the twist 3 correlator is related to causality property 
and compared with the prescriptions for exclusive hard processes. As a result we get the extra 
contributions, which naively do not have an imaginary phase. The single spin asymmetry for the 
i Drell-Yan process is accordingly enhanced by the factor of two. 
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I. INTRODUCTION 



Recently, the problem of the electromagnetic gauge invariance in the deeply virtual Compton scattering 
(DVCS) and similar exclusive processes has intensively been discussed in the literature, see for example [IHEj]. 
Qh) This development explored the similarity with the earlier studied inclusive spin-dependent processes @, and 
the transverse component of momentum transfer in DVCS corresponds to the transverse spin in DIS. 

The gauge invariance of relevant amplitudes is ensured by means of twist three contributions and the use 
of the equations of motion providing a possibility to exclude the three-particle (quark-gluon) correlators from 
the amplitude. After combining with the two-particle correlator contributions, one get the gauge invariant 
expression for the physical amplitude or, in the case of lepton-hadron processes, for the corresponding hadron 
tensor Q. 

This method was originally developed in the case of the particular inclusive processes with transverse polarized 
hadrons, like structure function g2 in DIS Q and Single Spin Asymmetry (SSA) Q due to soft quark (fermionic 
poles [8]). At the same time, the colour gauge invariance of the so-called gluonicpoles contributions [9| was 
previously explored by other methods relying on the Wilson exponentials [IT| - il4 |. 

Here we combine the approaches described above and apply them in the relevant case of the Drell-Yan (DY) 
process where one of hadrons is the transversally polarized nucleon. 

The SSA in the DY process was first considered in QCD in the case [H, EH of the longitudinally polarized 
hadron. This observable is especially interesting if the second hadron is a pion, because of the sensitivity l?], 



18] 
21] 



to the shape of pion distribution amplitude, being currently the object of major interest [19|, |20j (see also 
and Ref. therein). 

The imaginary phases in the SSA with longitudinally polarized nucleon are due to the hard perturbative gluon 
loops [H, [l6[ or twist 4 contribution of the pion distribution amplitude [13, [H, [H| • At the same time, the source 
of the imaginary part, when one calculates the single spin asymmetry associated with P+P'^ — > ££ + X process, 
is the quark propagator in the diagrams with quark-gluon (twist three) correlators. This leads [23j to the gluonic 
pole contribution to SSA. It has been reproduced (up to the derivative term, corresponding to the case of single 
inclusive Drell-Yan process, when only one of the leptons is observed) in the case of the non-zero boundary 
condition imposed on gluon fields, and the asymmetric boundary conditions have been considered as a privileged 
ones ITil . The reason is that these boundary conditions provide the purely real quark-gluon function B v (xi, 22) 
which parameterizes ( , ip^f + A^ip) matrix element. By this fact the diagrams with two-particle correlators do not 
contribute to the imaginary part of the hadron tensor related to the SSA. This property seems quite natural, 
as the corresponding diagram does not have a cut capable of producing the imaginary phase (25[. 

In our paper, we perform a thorough analysis of the transverse polarized DY hadron tensor in the light of 
the QED gauge invariance, the causality and gluonic pole contributions. 
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We show that to restore the electromagnetic gauge invariance of the transverse polarized DY hadron tensor, 
it is mandatory to add the extra diagram contribution (cf. [Ill), also at the twist three level. In contrast to 
the naive assumption, we demonstrate that this new additional contribution is directly related to the certain 
complex prescription in the gluonic pole 1 /(cci — x 2 ) of the quark-gluon function B (x±, x 2 ). It is essential that 
thisprescription is process-dependent, supporting the idea of effective process-dependent Sivers function (see 
e -g-[27[ and Refs. therein) related to this correlator. 

In more detail, we show that the causal pole prescription in the quark propagator, involved in the hard 
part of the standard diagram, supports the choice of a contour gauge and, in turn, the representation of the 
quark-gluon function B (xi, x 2 ) m the form of the gluonic pole with the mentioned complex prescription. This 
representation must be extended on the diagram, which naively does not contribute to the imaginary part. 
They ensure a new contribution to the imaginary part which is necessary to maintain the electromagnetic gauge 
invariance. Finally, the account for this new contributions corrects the SSA formula for the transverse polarized 
Drell-Yan process by the factor of 2. 



II. CAUSALITY AND CONTOUR GAUGE FOR THE GLUONIC POLE 



We study the contribution to the hadron tensor which is related to the single spin (left-right) asymmetry 
measured in the Drell-Yan process with the transversely polarized nucleon: N^^(p±) + N(p2) — >• 7*(<z) + 
X(Px) — > £{W) + £{h) + X(Px), where the virtual photon producing the lepton pair (li + 1 2 = q) has a large 
mass squared (q 2 = Q 2 ) while the transverse momenta are small and integrated out. The left-right asymmetry 
means that the transverse momenta of the leptons are correlated with the direction Sxe z where implies 
the transverse polarization vector of the nucleon while e z is a beam direction [28j . 

The DY process with the transversely polarized target manifests [23| the gluonic pole contributions. Since 
we perform our calculations within a collinear factorization, it is convenient (see, e.g., [29| ) to fix the dominant 
light-cone directions for the DY process shown at Fig. [T] 
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with n* = (1/V2, T , 1/V2), n M = (1/V2, T , -1/V2) 



(1) 



so that the hadron momenta p\ and p 2 have the plus and minus dominant light-cone components, respectively. 
Accordingly, the quark and gluon momenta ki and £ lie along the plus direction while the antiquark momentum 
k 2 - along the minus direction. 

Focusing on the Dirac vector projection, containing the gluonic pole, let us start with the standard hadron 
tensor generated by the diagram depicted on Fig. (Ha): 



WW 



d*k 1 d*k 2 5W(k 1 +k 2 -q) / d i £<S>^ )[ ^\k u £)^\k 2 )y. 



tr 



7m7 7*7 7< 



£ + j — k 2 7 + 



~2£+K 



where 



$Mh + l(ki,£) = 2 ( Pll S T \^( m h + gA a (z)^(0)\S T lPl ), $fr £ (paM»/2)7"^(0)ba 
Throughout this paper, J-\ and T 2 denote the Fourier transformation with the measures 

d i me ik2 - r > 2 and d% d*Z e - ik ^- u ' z , 
respectively, while and T 2 X mark the inverse Fourier transformation with the measures 

dye iyX and d Xl dx 2 e ixlXl+i(x2 ~ x ^ X2 . 



(2) 

(3) 
(4) 
(5) 



Analyzing the 7-structure of ((2|), we may conclude that the first term in the quark propagator singles out the 
combination: 7 + 7 Q 7~ with a = T which will lead to the matrix element of the twist three operator, (i/> ^A^tfj) 
with the transverse gluon field. After factorization, this matrix element will be parametrized via the function 
B (xi, x 2 ). The second term in the numerator of the quark propagator separates out the combination 7 + 7 a 7 + 
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with a = — . Therefore, this term will give (ip 7 + A + tp) which, as we will see now, will be exponentiated in the 
Wilson line [— oo~, CP]. Indeed, this part of the standard hadron tensor is given by 



w (l) [fe2 --tcnn] = , Mk .. Xl)y)tl 



-Id* 



7^7 7u7 7 7 

-il + z 



£+ -ie 



$ [7 ~ ] (fc 2 )x (6) 
d% e"*" 71 (pi, S T \i>( m ) 7 + gA+(0, z~ , 6 T ) m\S T ,Pi) , 



where 



dfj,(ki;xi,y) = dxid 4 kiS(xi - /pf) dyd 4 k 2 S(y - k 2 /p 2 ) 



S (4) (xipi + yp 2 - q) 



(7) 



Note that the prescription — ie in the denominator of this expression directly follows from the standard (see, 
e.g. [35j ) causal prescription for the massless quark propagator in 

Integrating over £ + , one can immediately obtain the corresponding ^-function in ©: 



yy (1) [ fc 2 -term] 



d/j,(ki;xi,y)tT 



7^7 7^7 



^ ] (fc 2 ) x 



+ 00 



e"^-" 1 (pi, S' i '|^(r 7 i)7 + ifl / dz~ 6(- Z -)A + {0, z~, T ) V(0)|5 T ,pi) 



(8) 



Including all gluon emissions from the antiquark going from the upper blob on Fig. HJa) (the so-called initial 

state interactions), we get the corresponding P-exponential in ^ ' (ki,£). The latter is now represented by 
the following matrix element: 

dVe- 4 * 1 -^!, S T m m ) 7 + [-cx)-, 0-]^(0)|5 T , Pl ), (9) 



where 



o 

[ -»-,0-] = Pexp{- <s /, i ,-^(0,,-,O T )}. 



(10) 



If we include in the consideration the gluon emission from the incoming antiquark (the mirror contributions) , 
we will obtain the Wilson line \rf[ ,— oo~] which will ultimately give us, together with (fT0|) . the Wilson line 
connecting the points and 771 in (|9|). This is exactly what happens, say, in the spin-averaged DY process [30l j . 
However, for the SSA, these two diagrams should be considered individually. Indeed, their contributions to 
SSAs, contrary to spin-averaged case, differ in sign and the dependence on the boundary point at — oo~ does 
not cancel. 

To eliminate the unphysical gluons from our consideration and use the factorization scheme @ , we may choose 
a contour gauge [3l| 

[-oo-0-] = l (11) 

which actually implies also the axial gauge A + = used in 

Let us discuss the problem of gauge choice in more detail. In pT|) . the so-called starting point xo (see, (3lj ]) 
is hxed to be at — oo~ owing to the certain complex prescription +ie in the quark propagator in ([2J. If we 
would change the starting point xq on +00 ~, this would correspond to the choice of the "anticausal" complex 
prescription —ie. On the other hand, the axial gauge A + = is independent on the choice of Xq and we are able 
to eliminate the Wilson line by choosing simply A + = without referring to the starting point x$. Nevertheless, 
since our prescription +ie in the quark propagator uniquely fixes the starting point Xq at —00, the expression 
for the Wilson line (fTTJj) hints the choice of gauge pT|) . 

Imposing this gauge one arrives [3l| at the following representation of the gluon field in terms of the strength 
tensor: 



(12) 
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Moreover, as we will demonstrate below, if we choose instead an alternative representation for the gluon in the 
form: 

oo 

A fi (z) = - { du)-6(u- - z-)G+^(u)-) + A"(oo) (13) 



(which corresponds to the gauge condition [+oo~, 0~] = 1 and also results in A + — 0) keeping the causal 
prescription +ie in ([2]), the cost of this will be the breaking of the electromagnetic gauge invariance for the DY 
tensor. 

We are now ready to pass to the term with £ + j~ in (j5J| which gives us finally the matrix element of the twist 
three quark-gluon operator with the transverse gluon field. Let us stop, in more detail, on the parametrization 
of the relevant matrix elements: 

( Pll S T \^\ 1 n) 1 pgA T a {\ 2 n)^)\S T , Pl ) ^ ie paS T pi B v (x u x 2 ) . (14) 
Using the representation (| 12[) . this function can be expressed as 

B v {x u x 2 ) = T ^ x *\ + S ( Xl - x 2 )B% ( Xl ) , (15) 

X\ — X 2 T tc 

where the real regular function T(x\, x 2 ) ( T{x, x) ^ 0) parametrizes the vector matrix element of the operator 
involving the tensor G M1/ (cf. [SH): 

{pi,S T \^(\ih)^ph v G va (\ 2 n)ip{ti)\S T ,pi) = e paS T pi T(xi,x 2 ) . (16) 
Owing to the time-reversal invariance, the function > (xi), 

izpaSTpt S(x 1 - x 2 )B\ (±oo) {x 1 ) = (pi, S T \ip(\in) -yp gA^{±oo) ip(0)\S T ,p{) , (17) 
can be chosen as 

Sa(-oo)(*)=0- (18) 

Indeed, the function B v (xi, x 2 ) is an antisymmetric function of its arguments @, while the anti-symmetrization 
of the additional term with ^X(-oo)( a; i) gi yes zero. 

There is no doubt that the only source of the imaginary part of the hadron tensor is the quark propagator. 
One may try to realize this property by assumption that matrix elements are purely real, 

B v (x u x 2 )= - T(xi,x 2 ), (19) 
x\ — x 2 

corresponding to asymmetric boundary condition for gluons (24|: 

B V A{oo) {x) = -B v A{ _ oo) {x) (20) 

Here we suggest another way of reasoning. The causal prescription for the quark propagator, generating its 
imaginary part, simultaneously leads to the imaginary part of the gluonic pole. Let us emphasize that this 
does not mean the appearance of imaginary part of matrix element (which by itself does not have an explicit 
physical meaning) but rather the prescription of its convolution with hard part. This procedure is in agreement 
with the prescriptions which were appeared in the exclusive case in the parametrization of the generalized gluon 
distributions |33l l34j|. 

Note that the fixed complex prescription +ie in the gluonic pole of B v (xi,x 2 ) (see, (fTS"]) ) is one of our main 
results and is very crucial for a new contribution to hadron tensor we are now ready to explore. Indeed, the 
gauge condition must be the same for all the diagrams, and it leads to the appearance of imaginary phase of 
the diagram (see, Fig. [ljb)) which naively does not have it. Let us confirm this by explicit calculation. 
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Figure 1: The Feynman diagrams which contribute to the polarized Drell-Yan hadron tensor. 



III. HADRON TENSOR AND GAUGE INVARIANCE 

We now return to the hadron tensor and calculate the part involving £ + ^/~ , obtaining the following expression 
for the standard hadron tensor (see, the diagram on Fig. da)): 



— (l)[*+-tcrm] 



/ dx2 tr 



= J d 2 <i T W$ = -J d Xl dy 

T (%1 ~ X 2 )p\ 



S(xi - x B )8(y - y B ) 



q(y) x 



(21) 



[xi - x 2 )ys + ie 



B V (xi,X 2 ) £p a S T Pl 7 



where we used £ +r ) = (x 2 — xi)pi and 



(p 2 1 1/>( An*) 7^^(0)1^2) = p 2 ^q{y)- 



(22) 



We are now in position to check the QED gauge invariance by contraction with the photon momentum q^. 
Calculating the trace 



1 



(xi - x 2 )ep aS T pi tr 



qjpl„p 2 l a Pi 



(23) 



one gets 



In w£2 = - I dxi dy 



5{x\ - x B )S(y - y B ) 



q(y)£v P2 ST pi I dx 2 3m— — ^—B v (x 1 ,x 2 )^0, (24) 

xi — x 2 + it 



if the gluonic pole is present. Note that here and below we consider only the imaginary part of the hadron 
tensor ( as for any single spin asymmetry). 

Let us analyze this problem from a viewpoint of the so-called ^-process (see (35[, Section 33.2) applied for 
the partonic sub-process. Generally speaking, the single diagram on Fig. HJa) cannot give the gauge invariant 
hadron tensor. One needs the second diagram (cf. ) with the gluon insertion in the quark line, see Fig. [Ub). 

We now focus on the contribution from the diagram depicted on Fig. []Jb). The corresponding hadron tensor 
takes the form: 



7 A1 J r (fci)7 1/ $(/c 2 ) 



W$ = f d 4 ^ d 4 k 2 5^{h + k 2 - g)tr 
where the function .F(fei) reads 

F{kx) = Sfaft* I dSie- lk ^( Pl ,S T \i ! ( m )gAl(0)i;(0)\S T ,Pi) 



(25) 



(26) 



Performing the collinear factorization, we derive the expression for the factorized hadron tensor which corre- 
sponds to the diagram on Fig. Hfb): 



w (2) 



dx\ dy 



S(xi - x B )5(y - y B ) 



q{y) tr 



7 M ( / d 4 fci S(xipi - ki)F(ki) I71/P2 



(27) 



G 



After some algebra, the integral over k\ in (j2"T)l can be rewritten as 



j d A k x 5(x lP + - kt)T^ + \k x ) = ^3£^ e( ,. s , 



2p 2 p{ xx + ie 



dX2 B V (x 1 ,X2) ■ 



(28) 



where the parametrization (fT4|) has been used. Taking into account (|28|) and calculating the Dirac trace, the 
contraction of the tensor Wp V with the photon momentum gives us 



9m = da;i ^ 



S(xi - x B )5{y - Vb) 



1 



(29) 



From this, one can observe that if the function B v (x 1,2:2) is the purely real one (see, (fl9| ). this part of the 
hadron tensor, which is associated with the diagram on Fig. [3Jb), does not contribute to the imaginary part. 

We now study the net effect of the VvA and wj^ contributions and its role for the QED gauge invariance. 
adding the contributions of (|24l) and (|29| . one can easily obtain: 



9m W m^ + 9m 



~-vp2S T ]: 



pl q(yB)^m dx 2 B (x B ,x 2 ) 



x B - x 2 
xb - x 2 + 



- 1 



(30) 



If we tacitly assume that B v (xi, x 2 ) is some real and regular (at x\ = x 2 ) function that the numerator and 
denominator in the first term inside the brackets are contracted and, as a result of this, both the first and 
second terms in (|30|) do not have an imaginary part. That would mean the electromagnetic gauge invariance 
for the tensor. 

The existence of the gluonic pole changes the situation. Inserting now (fT5|) into (|30"|) , one gets 



9m W2 



9m VlC 



£ vP2S T Pl 9(ys)^m / dx 2 T(x Bl x 2 ) 
-1 



SB - X 2 



1 



(xb - x 2 + ie) 2 



xb - x 2 



(31) 



Performing the calculation one gets: 



9m 



(2) 



0. 



(32) 



This is nothing else than the QED gauge invariance for the imaginary part of the hadron tensor. From (|31|) . 
we can see that the gauge invariance takes place only if the prescriptions in the gluonic pole and in the quark 
propagator of the hard part are coinciding. Indeed Eq. pip with the field (fT3")) takes the form 

9m + g„ = 



£ vp2 S T P i q{Vb) 3m / dx 2 T(x B ,x 2 ) 
-1 



x B - x 2 



1 



(a;_B — x 2 — ie)(xB — x 2 + ie) xb — x 2 + ie 



(33) 



It is clear that the first term in the brackets is purely real, and the imaginary part from the second term stays 
uncompensated. Let us note for completeness, that the treatment of the pole in the principal value sense is 
equivalent to the mean arithmetic of two discussed prescriptions and also cannot satisfy the gauge invariance. 
Thus we completed the reductio ad absurdum of the hint suggested in Section [TT] and found that the contour 
gauge (jlip is a correct one. In other words, it means that the prescription in the quark propagator must agree 
with the representation of B v (xb, x 2 ). Otherwise, one may face the problem with the gauge invariance. 

It is instructive to compare the electromagnetic gauge invariance of the gluonic poles contributions with that 
of perturbative QCD. In the latter case the imaginary part is provided by hard gluon loops and the QED gauge 
invariant set consists of 3 diagrams depicted on Fig. [2] At the same time, the imaginary part is due to the 
single diagram on Fig. [2ja) and it is gauge invariant by itself as the photon line couples to two on-shell (because 
of the Cutkosky cutting rule) quarks. This reasoning, however, does not happen to work for (non-perturbative) 
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Figure 2: The Feynman diagrams which contribute to the a -order amplitude in QED. 



gluonic pole contribution (see, Fig. Ufa)) an d the contribution of the diagram on Fig. [TJb) should be added 
to ensure the electromagnetic gauge invariance. This is clearly seen from Eq. ([30]) where the analog of the 
contribution of the diagram on Fig. Eta) is represented by the first term in the brackets. Its imaginary part 
is zero (i.e. QED gauge invariant) only if gluonic pole is absent at all. This situation corresponds also to the 
difficulties in the applicability of Ward identities to gluonic poles contributions (see [36| and Refs. therein). 

As we have shown, only the sum of two contributions represented by the diagrams on Fig. QJa) and (b) can 
ensure the electromagnetic gauge invariance. We now inspect the influence of a "new" contribution QJb) on the 
single spin asymmetry and obtain the QED gauge invariant expression for the hadron tensor. It reads 



yrfll = wj^ + y^l = -^2 £ »s T P1P2 z v i(vb) t(x b ,x b ) 



where one used q 2 = sxbUb and introduced the vector 



Zp, = pi M - P2 fj, = XB Pi n - UB P2 f 



which together with the vectors: 



Y, 



form the mutually orthogonal basis (see, [28J). Here pi M are the partonic momenta (q^ =pi^ 
help of ([33]) and (j3"6")l . the lepton momenta can be written as (this is the lepton cm. system) 



where A = A/ \/— A 2 and 

/ M (#, <p; X, Y, Z) — coscp sin# + F M sin<^ sin0 + Z^ cosi 



Z) 



(34) 
(35) 
(36) 

-p 2fl ) With the 
(37) 

(38) 
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Within this frame, the contraction of the lepton tensor with the gauge invariant hadron tensor reads 

QJ 

^VV^ = -2cos6»e iyS T pip2 q(y B )T(x Eh x B ) • (39) 

We want to emphasize that this differs by the factor of 2 in comparison with the case where only one diagram, 
presented on Fig. HJa), has been included in the (gauge non- invariant) hadron tensor, i.e. 

c»uW { tl = \^»wfu- (40) 

Therefore, from the practical point of view, the neglecting of the diagram on Fig. [TJb) or, in other words, the 
use of the QED gauge non-invariant hadron tensor yields the error of the factor of two. 

Indeed, taking the contribution of the diagram Fig. QJa) corresponds to keeping of only the term proportional 
to pi fj, in (|35|) . The contraction with (gauge invariant) lcptonic tensor is equivalent to making it gauge invariant 
by substitution 



P_q = PiMjvPjjf 
1 Q 2 2 

It is this factor of 2 which makes the difference with the correct gauge invariant expression. 



Pi/* =>Pi/* ~ %~Q2 = 2 ' ^ ' 



IV. CONCLUSIONS AND DISCUSSIONS 



The essence of this paper consists in the exploration of the electromagnetic gauge invariance of the transverse 
polarized DY hadron tensor. We showed that it is mandatory to include a new contribution of the extra 
diagram which naively does not have an imaginary part. The account for this extra contribution leads to the 
amplification of SSA by the factor of 2. 

This new additional contribution emanates from the complex gluonic pole prescription in the representation 
of the twist 3 correlator B {x\, x%) which, in its turn, is directly related to the complex pole prescription in 
the quark propagator forming the hard part of the corresponding hadron tensor. 

In more detail, the causal prescription in the quark propagator, involved in the hard part of the diagram on 
Fig. [Ha), selects from the physical axial gauges the contour gauge defined by Eq. (fTTj) . At the same time, the 
contour gauge predestines Eq. (|12p and, therefore, the representation of B {x\,X2) in the form of the gluonic 
pole with the complex prescription, see (|15[) . Since both diagrams on Fig. QJa) and (b) should be considered 
within the same (contour) gauge, the representation (|15[) . which we advocate, has to be applied for the diagram 
depicted on Fig. [TJb). As a result of this, the diagram on Fig. [ljb), in contrast to naive assumptions, has the 
imaginary part. In some sense, the diagram on Fig. [IJb) feels the complex prescription in the hard part of the 
diagram on Fig. HJa) by means of the contour gauge which we make used. Note that, from the physical point 
of view, the consideration of each of the diagrams on Fig. Q] individually makes no sense. 

This is completely similar to the case of exclusive dijet production [34j when the pole prescription in (twist 
two) matrix element of gluonic fields is controlled by the corresponding hard subprocess. 

We have argued that, in addition to the electromagnetic gauge invariance, the inclusion of new-found contri- 
butions corrects by the factor of 2 the expression for SSA in the transverse polarized Drell-Yan process. 

Finally, we proved that the complex prescription in the quark propagator forming the hard part of the hadron 
tensor, the starting point in the contour gauge, the representation of B (xi, x<z) like tjl 5|) and the electromagnetic 
gauge invariance of the hadron tensor must be considered together as the deeply related items. 
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